February 22, 2008 14:15 



WSPC - Proceedings Trim Size: llin x 8.5in reville 



Environmental limits on the non-resonant cosmic-ray current-driven instability 



B. Reville, J.G. Kirk 
Max- Planck- Institut fur Kernphysik, Heidelberg 69029 



oo 

o 

O 

(N 

X> 

(N 
(N 



43 

Or 
i 

o 



> 

(N 

m 
m 

(N 
O 
OO 

o 



X 



P. Duffy 

UCD School of Physics, University College Dublin, Belfield, Dublin 4 
S. O'Sullivan 

School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9 



We investigate the so-called non-resonant cosmic-ray streaming instability, first discussed by Bell (2004). The extent 
to which thermal damping and ion-neutral collisions reduce the growth of this instability is calculated. Limits on the 
growth of the non-resonant mode in SN1006 and RX J1713. 7-3946 are presented. 
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1. Introduction 

Diffusive shock acceleration at the outer shocks of 
supernova remnants is believed to be the primary 
source for galactic cosmic rays. In the simplest pic- 
ture, however, the maximum energy attainable by 
this mechanism falls short of the so-called knee at 
approx. 1 PeV [1]. Amplification of the magnetic 
field in the vicinity of shock front may facilitate ac- 
celeration to higher energies, without changing the 
simplicity of the model. There has also been an in- 
crease in the wealth of observational evidence for 
magnetic field amplification at collisionless shocks, 
such as the non-thermal bright X-ray rims in several 
supernova remnants [2, 3]. The growth of the mag- 
netic field is thought to be closely related to the effi- 
cient production of cosmic rays and the non-resonant 
current-driven instability [4] has provided a mech- 
anism which may indeed considerably amplify the 
magnetic field. 

Generally speaking, the medium upstream of a 
supernova blastwave will not be completely ionised. 
The presence of a small neutral component limits the 
acceleration of cosmic rays in these environments [5] , 
due to damping of resonantly excited Alfven waves. 
In molecular clouds, where the ionisation fractions 
are even lower and the densities much higher, the 
situation is even more extreme [6, 7]. In a certain 
range of wavelengths, Alfven waves simply do not 
propagate, and are quite rapidly damped. 

A non-resonant mode, driven by a cosmic-ray 
current in a partially ionised medium has recently 



been studied [8]. Treating the background particles 
as a MHD single fluid with a generalised Ohm's law, 
and using a kinetic description of the cosmic rays, 
these authors conclude that a non-resonant mode 
with a rapid growth rate remains. However, as we 
show, for high frequency waves the friction between 
the charged species and the neutrals is not sufficient 
to couple their motions, and a single fluid description 
of the plasma is no longer appropriate. 

In the present work, we investigate various condi- 
tions that may reduce the growth of the non-resonant 
instability. We first calculate the effect of thermal 
damping of the waves. We also determine an ex- 
pression for the dispersion relation of waves in the 
presence of an arbitrary fraction of neutral particles. 
These results are then used to speculate on whether 
or not magnetic field amplification may have taken 
place in SN1006 or in the molecular cloud close to 
RX J 1713. 7-3946. We conclude with some remarks 
about the importance of the non-resonant current- 
driven instability. 

2. Thermal effects 

We present here a simple derivation of the cosmic-ray 
current-driven instability. Since the fastest growing 
mode is non-resonant we do not need to calculate the 
full effects of an anisotropic cosmic-ray distribution 
[9], and we investigate here a super- Alfenic beam, 
which considerably simplifies the analysis. 

The linear dispersion relation for circularly po- 
larised transverse waves propagating parallel to the 
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zeroth order magnetic field is 
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where the summation is over each species with cor- 
responding charge q s , cyclotron frequency w cs = 
q s B /m s c and plasma frequency lo 2 s — 4Trq 2 n s /m s . 
The susceptibility Xs for each component of a plasma 
is determined by integrating the Vlasov equation 
along the unperturbed trajectories about the zeroth 
order field and can be found in any standard text- 
book on plasma physics (e.g. Krall & Trivelpiece 
1973). For circularly polarised waves propagating 
parallel to a uniform magnetic field, the susceptibil- 
ity for each species is given by 
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For simplicity, the plasma we consider has only 
three components. A two species thermal population 
of electrons and protons, and a cosmic ray compo- 
nent which we assume to consist entirely of pro- 
tons. Since we are investigating waves on length- 
scales longer than the Debye length the plasma can 
be taken to be quasi-neutral, resulting in a slight 
charge excess of electrons in the thermal plasma due 
to the presence of cosmic rays. A neutralising back- 
ground flux is also induced to balance the current 
produced by the streaming cosmic rays. 

Assuming a Maxwellian background and a 
mono-energetic beam of protons with momentum po 
directed along the magnetic field, it is straightfor- 
ward to show that the dispersion relation for low 
frequency waves |o»| <C uJ c i < |w ce | is 
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where c/3b is the speed of the proton beam, va = 
Bo/^Annimi is the Alfvcn velocity and 6 = 
k^Ti/miC 2 the dimensionless ion temperature. We 
have also introduced the parameter e to describe the 
polarisation with e = +1(— 1) for right (left)-handed 
waves (cj > 0). 

We define r g = pqc/cBq the gyroradius of the 
protons in the beam and consider only waves with 
k > 0. We have also introduced a dimensionless pa- 
rameter characterising the strength of the driving 
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The maximum growth rate of the non-resonant 



mode is 



Im(w) w (3 h u c 

2 va rii 



(4) 



which is independent of magnetic field strength. 

We can now determine possible thermal effects 
on the growth rate. For kr g 3> 1, it follows from (3), 
in the limit 8>va 2 /c 2 , that the maximum growth 
rate is 



(5) 



In order for the non-resonant mode to leave the 
regime of linear growth before being overtaken by 
the shock front, i.e., before being advected over a dis- 
tance of roughly r g //3b, one requires Im(o;) > (3 2 c/r g . 
The necessary condition for thermal effects to reduce 
the growth rate below this value is 

C 2 Po 



e > 



(6) 



For typical SNR parameters this condition will only 
be satisfied for a very weak driving term , ( < 1. 
However, thermal effects are likely to play a signif- 
icant role for relativistic shocks, and may even pro- 
vide a saturation mechanism for the current driven 
instability [10]. 

3. Collisional effects 

The surroundings into which supernova shocks prop- 
agate are, in general, not completely ionised. Obser- 
vations have shown that many supernova take place 
near molecular clouds, which can have very low ion- 
isation fraction. We consider the reduction of the 
growth of the non-resonant modes due to ion-neutral 
friction. For wave frequencies larger than the momen- 
tum exchange frequency, a multi-fluid treatment of 
the plasma is necessary, since the collisions are not 
frequent enough to couple the motion of the neutral 
particles to the charged ions. We consider two fluids, 
a single MHD fluid for the charged species and a neu- 
tral component. A simple analysis of the frictionally 
coupled system with an external cosmic ray current 
results in the following dispersion relation [9] 

,,2 



1 + 



= k 2 v A 2 + e(^k(a-l). (7) 



February 22, 2008 14:15 WSPC - Proceedings Trim Size: llin x 8.5in reville 



3 



where v a b is the momentum exchange frequency from 
species a to species b and are related by momen- 
tum conservation p n v n i — p\y m . The complex func- 
tion cr(fc) represents the susceptibility of the cosmic 
rays normalised to the induced thermal return cur- 
rent [4, 10]. For plasmas in the temperature range 
10 2 K < T < 10 5 K, [11] give the following expres- 
sion 



8.9 x l(T 9 n n 



T 



10 4 K 



In the high frequency limit, \u>\ 3> v\m this re- 
duces to a dispersion relation similar to Eq. (3). In 
the low frequency limit \w\ <C iy nl we find 



u) i=a — 



A: 2 v A 2 + eC— k(a-l) 



(9) 



If the collision frequencies are very large, the neutral 
and ionised components are tied together and the 
effect of ion-neutral collisions is simply to increase 
the effective mass of the ions. 

For strongly driven, non-resonant modes Eq. (7) 
reduces to 



uj(oj 2 +(— k)+iv n 

Tg; 



(1 + Z)u? + Z(^k 



0, (10) 



where Z = pi/p n = fni/^in- For all physically rel- 
evant parameters, this cubic polynomial in w has 
three purely imaginary roots, two of which are 
damped modes, and the third the non-resonant grow- 
ing mode. Thus, ion-neutral collisions are unable to 
stabilise the strongly driven mode, although they can 
affect its growth rate. Interestingly this result also 
applies to any aperiodic mode. 

The general analytic expressions for the grow- 
ing mode are cumbersome, but various limiting cases 
yield interesting results. In the limit of low ionisation, 
Z ^> 1 one recovers the growth rate given in Eq. (4). 
For Z <C 1 the growth rate is given by 



7 = — 




(11) 



If v? n <C 4Cv 2 fc/r g collisions are too slow to compete 
with the driving, and the growth rate is again the 
same as Eq. (4) with a negligible reduction due to 
collisions. However, in addition to the low frequency 
limit in which the components are tied together Eq. 
(9), a new regime arises. If the collisional drag on 
the ions is sufficient to affect the driving, but not to 



couple the components, i.e., uf n ^> 4£v 2 /c/r g 3> v 2 -^ 



then we obtain 



7«C- 



■fcv s . 



(12) 



'in ' g ^in Hi 

A similar result to this has been found previously [8] , 
but as we have shown here, this is valid only when 

7 < ^in- 



(8) ImH 
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Fig. 1. Maximum growth rate as a function of ionisation 
fraction for total density n = n\ + n n in a H-H + gas: 
0.1 cm- 3 , 10 4 K (solid), 1.0 cm" 3 , 10 3 K (dash), 10.0 cm" 3 , 
10 2 K (dash-dot). We take £ = 0.01xiv s /c with shock speed 
v s = 5000 km s — 1 , and r g = 10 14 cm. ui is in units of v 2 /r g c. 
Once the o; max < 1 the non-resonant instability is not effective 
in amplifying the magnetic field. 



In Fig. 1 we illustrate the influence of ion-neutral 
collisions for SNR parameters using a full numerical 
solution to the dispersion relation. Thermal damping 
is negligible, but a weak dependence on temperature 
enters via Eq. (8). We plot the maximum growth 
rate, as a function of ionisation fraction X[ = rii/n, 
taking the driving term £ = 0.01x;v s /c since the cos- 
mic rays couple directly only to the ionised compo- 
nent. For higher density plasmas, a larger fraction of 
neutral particles is necessary to reduce the growth 
rate below the threshold value of v 2 /r g c. 

4. Application to SNR 

We now investigate the role of the non-resonant 
instability in some typical supernova remants. Al- 
though it is beyond the scope of this work to esti- 
mate the possible saturated field strength from the 
non-resonant instability, we can nevertheless test the 
consistency of the analysis, by checking that the in- 
stability can enter the non-linear regime before being 
overtaken by the shock. Thus, the physical conditions 
far upstream (£?, n, x- u T) together with the shock 



February 22, 2008 14:15 



WSPC - Proceedings Trim Size: llin x 8.5in reville 



speed v s and cosmic ray intensity U CY /piV s 2 = 0.1, 
should combine to yield a growth rate in excess of 
v s 2 /(r g c). We investigate two extreme cases. 

Firstly, for SN1006, the shock velocity is approx- 
imately v s w 2900 km s~ the density is 0.05 < n < 
0.25 cm -3 , and the neutral fraction is small, x\ w 0.9 
[12]. Taking m = 0.25 cm" 3 gives Cv s 2 / v a 2 ~ 10. 
The low density and shock speed, imply that the non- 
resonant mode is not very strongly driven, but the 
driving is nevertheless strong compared to the colli- 
sion frequency. Very close to the shock, the growth 
time of the field is on the order of years, while fur- 
ther from the shock the growth rate decreases quite 
dramatically, as the number density of cosmic rays 
decreases. For the parameters given above the max- 
imum growth rate is ~ 10v s 2 /r g c. The non-resonant 
mode may still amplify the field into the nonlinear 
regime. 

The low neutral fraction in SN1006 means that 
ion-neutral collisions are not likely to dramatically 
reduce the instability. The interaction of a supernova 
blast wave with a molecular cloud, such as that ob- 
served close to the north-western rim of RX J1713.7- 
3946 would have a much more significant effect on 
the growth. Although the exact parameters are un- 
certain, molecular clouds are generally clumpy with 
interclump densities in the range 5 — 25 cm~ 3 and 
ionisation fractions not larger than 10% [13]. Adopt- 
ing the values in [14] T = 10 2 K, n = 23 cm~ 3 , x { = 
0.01, v s = 10 8 cm s _1 the collision frequency is larger 
than the growth rate. The non-resonant driving term 
is unable to dominate over the Alfven term and 
growth is driven by the cosmic rays themselves, at 
very long wavelengths kr g <C 1. At shorter length 
scales, the waves are rapidly ion-neutral damped 
Alfven waves. The shortest growth timescale is on 
the order of ~ 10 6 yrs. This value varies with choice 
of r g , temperature and density, but for high density 
and low ionisation fractions, the growth timescale is 
longer than the free expansion phase of a supernova 
remnant, or the lifetime of the wave in the shock pre- 
cursor, suggesting that the Fermi acceleration mech- 
anism may switch off upon interaction with a molec- 
ular cloud. 

5. Discussion 

It is now widely believed that magnetic field amplifi- 
cation is the best method for pushing the maximum 



energy beyond the Lagage-Cesarsky limit [1]. The 
non-resonant instability has been the most promis- 
ing mechanism so far in that its growth rate is con- 
siderably faster than its resonant counterpart, and 
initial studies suggest that it does not saturate at 
a level SB < Bo. However, since the mechanism 
is non-resonant, with the fastest growing modes at 
very short wavelengths, the amplified field has very 
little effect on the acceleration of the maximum en- 
ergy particles unless the energy is transferred to the 
large scale field by some non-linear process, as seen 
in some numerical simulations [4, 15]. The results re- 
ported on in this paper extend the previous studies 
and allow application to more extreme environments. 
Such processes may be crucial for understanding par- 
ticle acceleration at supernova remnant shocks in the 
galactic centre region, for example. Efficient acceler- 
ation is essential for this mechanism to occur. The 
largest uncertainty in our model arises in determin- 
ing the value of the driving term £. It is suggested, 
however, that efficient acceleration naturally leads to 
magnetic field amplification. 
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